Critical behavior of the two-dimensional TV-component Landau-Ginzburg Hamiltonian 

with cubic anisotropy. 
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. inf n. it 



We study the two-dimensional TV-component Landau-Ginzburg Hamiltonian with cubic anisotropy. 
We compute and analyze the fixed-dimension perturbative expansion of the renormalization-group 
functions to four loops. The relations of these models with TV-color Ashkin- Teller models, discrete 
cubic models, planar model with fourth order anisotropy, and structural phase transition in adsorbed 
monolayers are discussed. Our results for N — 2 (XY model with cubic anisotropy) are compatible 
with the existence of a line of fixed points joining the Ising and the 0(2) fixed points. Along this line 
the exponent r\ has the constant value 1/4, while the exponent v runs in a continuous and monotonic 
way from 1 to oo (from Ising to 0(2)). For N > 3 we find a cubic fixed point in the region u, v > 0, 
which is marginally stable or unstable according to the sign of the perturbation. For the physical 
relevant case of TV = 3 we find the exponents 77 = 0.17(8) and v = 1.3(3) at the cubic transition. 



PACS Numbers: 75.10.Hk, 05.70.Jk, 64.60. Fr, ll.10.Kk 



I. INTRODUCTION 

In the framework of Renormalization Group (RG) ap- 
proach to critical phenomena, the critical behavior at 
many continuous phase transitions can be investigated by 
considering an effective Landau-Ginzburg Hamiltonian, 
having an TV-component order parameter fa as funda- 
mental field, and containing up to fourth-order powers 
of the field components. The fourth-degree polynomial 
form of the potential depends essentially on the symme- 
try of the system. In fact according to the universality 
hypothesis, the critical properties of these systems can be 
described in terms of quantities that do not depend on 
the microscopic details of the system, but only on global 
properties such as the dimensionality and the symmetry 
of the order parameter, and the range of the interactions. 

The critical properties of many magnetic materials are 
computed using the 0(TV)-invariant Landau-Ginzburg 
Hamiltonian. Uniaxial ferromagnets should be described 
by the Ising universality class (TV = 1), while magnets 
with easy-plane anisotropy should belong to the XY 
universality class. Ferromagnets are often described in 
terms of the 0(3) Hamiltonian. However, this is correct 
if the non-rotationally invariant interactions that have 
only the reduced symmetry of the lattice are irrelevant 
in the renormalization-group sense. In two dimensions 
the effect of anisotropy is very important: systems pos- 
sessing continuous symmetry do not exhibit conventional 
long-range order at finite temperature, while models with 
discrete symmetry do undergo phase transitions into con- 
ventionally ordered phase. 

For studying the effect of cubic anisotropies one usually 
consider the </> 4 theory 
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in which the added cubic term breaks explicitly the O(N) 
invariance leaving a residual discrete cubic symmetry 
given by the reflections and permutations of the field 
components (3|. This term favors the spin orientations 
towards the faces or the corners of an TV-dimensional hy- 
percube for vq < or vq > respectively. 

The Hamiltonian ||) has received much attention also 
because it describes the two-dimensional order-disorder 
transition in some adsorbed monolayers. In fact, in the 
original classification 0] on square and rectangular lat- 
tices, all these transitions belong to the universality class 
of the Ising model, three- or four-state Potts model, and 
XY model with cubic anisotropy. In the successive exten- 
sion ||] to non-Bravais lattices, as the honeycomb and the 
Kagome ones, also the Heiscnberg (TV = 3) model with 
cubic anisotropy is interested by this classification (e.g. 
in the adsorption of diatomic- molecules on graphite). 
The full classification of continuous phase transition of 
magnetic symmetry in two dimensions || reveals the in- 
terest of the model (Q) also for other kind of transitions. 

An important application of the XY model with cu- 
bic anisotropy is in the oxygen ordering in YBaCusOe+a; 

, since it is one of the most studied high-temperature 
superconductors. Some other applications of Landau- 
Ginzburg Hamiltonian (]l|) are: the buckling instabilities 
of a confined colloid crystal layer |J , some discrete mod- 
els with competing nearest- and next-nearest-neighbor 
interactions |hJ, and, very recently 11 1, it is argued that 
the phase diagram of a lattice-gas model for studying the 
micellar binary solution of water and ampliphile is par- 
tially described by the XY model with cubic anisotropy. 

We also mention that, in the limit TV — > 0, the cubic 
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model (|l|) describes the Ising model with site-diluted dis- 
order |12|-|l4|, but we do not consider it here, since in 
the framework of fixed dimension d = 2 it was already 
analyzed in Ref. | JL5| . 

In this paper we present a field-theoretic study based 
on an expansion performed directly in two dimensions, 
as proposed for the O(N) models by Parisi |16|. 

The paper is organized as follows. In Sec. II we give an 
overview of all known results that we believe necessary for 
a full understanding of the two-dimensional cubic model 
(|l|). In Sec. Ill we derive the perturbative series for the 
renormalization-group functions at four loops and discuss 
the singularities of the Borel transform. The results of 
the analysis are presented in Sec. IV. The reader that is 
not interested to the detail of the calculations can skip 
sections III and IV and read directly Sec. V where we 
summarize all our results and point out some questions 
which we think deserve further study. 



II. OVERVIEW OF KNOWN RESULTS 
A. Three dimensional results 

We shortly review the most interesting features ap- 
pearing in the study of the Hamiltonian (Q) in the frame- 
work of e-expansion []l7],[l| and at fixed dimension d = 3 
fl8|| , since they are necessary for a good understanding 
of the two-dimensional case. 

The model (|l|) has four fixed points: the trivial Gaus- 
sian one, the Ising one in which the N components of the 
field decouple, the 0(iV)-symmetric and the cubic fixed 
points. The Gaussian fixed point is always unstable, and 
so is the Ising fixed point p9[ . Indeed, in the latter case, 
it is natural to interpret Eq. ([[]) as the Hamiltonian of N 
Ising-like systems coupled by the 0(iV)-symmetric term. 
But this interaction is the sum of the products of the en- 
ergy operators of the different Ising systems. Therefore, 
at the Ising fixed point, the crossover exponent associated 
to the 0(iV)-symmetric quartic term should be given by 
the specific-heat critical exponent ai of the Ising model, 
independently of N. Since ai is positive for all d > 2 the 
Ising fixed point is unstable. Obviously in two dimen- 
sions this argument fails since ai = 0. 

While the Gaussian and the Ising fixed points are un- 
stable for any number of components N, the stability 
properties of the 0(7V)-symmetric and of the cubic fixed 
points depend on N. For sufficiently small values of N, 
N < N c , the 0(iV)-symmetric fixed point is stable and 
the cubic one is unstable. For N > N c , the opposite is 
true: the renormalization-group flow is driven towards 
the cubic fixed point, which now describes the generic 
critical behavior of the system. The 0(iV)-symmetric 
point corresponds to a tricritical transition. Figure |l| 
sketches the flow diagram in the two cases N < N c and 
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FIG. 1. Renormalization-group flow in the coupling plane 
(u, v) for N < N c and N > N c . 



N > N c . At N — N c , the two fixed points should coin- 
cide, and logarithmic corrections to the 0(iV)-symmetric 
critical exponents are expected. Outside the attraction 
domain of the fixed points, the flow goes away towards 
more negative values of u and/or v and finally reaches 
the region where the quartic interaction no longer satis- 
fies the stability condition. These trajectories should be 
related to first-order phase transition s f2(i[| . Some recent 
and very accurate calculations ]2l| p7| , ^8| suggest that in 
three-dimensions N c < 3. 

If N > N c , the cubic anisotropy is relevant and there- 
fore the critical behavior of the system is not described 
by the Heisenberg isotropic Hamiltonian. If the cubic 
interaction favors the alignment of the spins along the 
diagonals of the cube, i.e. for a positive coupling vq, the 
critical behavior is controlled by the cubic fixed point 
and the cubic symmetry is retained even at the critical 
point. On the other hand, if the system tends to magne- 
tize along the cubic axes — this corresponds to a negative 
coupling i>o — then the system undergoes a first-order 
phase transition ]j,g^-go| . 

In the limit N — > oo, keeping Nu and v fixed, one 
can derive exact expressions for the exponents at the cu- 
bic fixed point for all dimensions. Indeed, in this limit 
the model can be reinterpreted as a constrained Ising 
model jn]], leading to a Fisher renormalization of the 
Ising critical exponents J3^] . One has |3^|n],[l| : 



r\ = rji + 



1 — Oil 
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(2) 



where r)i, j//, and ai are the critical exponents of the 
Ising model. 

In all dimensions, a simple argument based on the sym- 
metry of the two-component cubic model |34| shows that 
the cubic fixed point for TV = 2 has the same stability 
properties of the Ising one. Indeed, for N = 2, a 7r/4 
internal rotation, i.e. 

{</>!> <fa) — ►-^=(01 + 02,01-02), (3) 
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maps the cubic Hamiltonian ([!]) into a new one of the 
same form but with new couplings (u' , v ) given by 



u' = u Q + |^o, 



-v . 



(4) 



This symmetry maps the Ising fixed point onto the cubic 
one. So for all d > 2, since the Ising point is unstable, 
the cubic point is unstable too, and the stable point is 
the isotropic one. In two dimensions, this is no longer 
true. Indeed, one expects the cubic interaction to be 
truly marginal for N = 2 HQ and relevant for N > 2 
[B7|, so N r = 2 in two dimensions. 



B. The two dimensional case 
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when hi — > 0. For this reason it is often said that the 
XY model with cubic anisotropy has a non-universal be- 
havior, although there is the so-called weak-universality 
of Suzuki |40| (i.e. r\ is constant). 

In a successive work Kadanoff EJ identifies the two 
fixed-point lines for /14 ^ with the dual line of 
the Ashkin- Teller (AT) and the eight-vertex Baxter 
model (8V) (see the Baxter book [EJI for a review about 
these models). One remark is necessary at this point: 
the AT and the 8V models are equivalent on the dual 
line E2) and they show a continuously varying exponent 
of the correlation length: 



As it is already clear from the previous subsection, in 
two-dimensions several new and interesting features ap- 
pear. First of all we have no general argument to under- 
stand the stability properties of the Ising fixed point, in 
fact the specific heat of the two-dimensional Ising model 
has a logarithmic divergence (i.e. ui — 0). The van- 
ishing of this crossover exponents, and so the presence 
of a marginal operator, can bring us to several different 
scenarios. For example we may have that the Ising fixed 
point is marginally stable or unstable because of higher 
order corrections to RG equations, or that there exists 
a line of fixed points (see Ref. 1 38 1 for a detailed review 
about the effects of marginal operators). 

When the cubic anisotropy becomes marginal (at N = 
N c = 2) |}7| the cubic fixed point is in the region with 
V < (cf. Eq. (Eh), and so it cannot coincide with the 
O(N) one, as in the case of three and 4 — e dimensions. 

There are several studies on systems related to Landau- 
Ginzburg Hamiltonian (fjl). In Ref. pa the model 
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was considered, where the first sum is only over nearest- 
neighbors and 9(r) is the angle that the two-component 
spin S(r) forms with some arbitrary axis. This sys- 
tem clearly belongs to the same universality class of XY 
model with cubic anisotropy if p = 4 (the "field" /14 maps 
on v Q ). 

The phase diagram found in Ref. ||^] consists of three 
distinct lines of fixed points showing continuously vary- 
ing exponents. The line with /14 = is the standard 
Kosterlitz-Thouless (KT) |§. There are two lines with 
hi =/= 0, starting from the KT transition (the end point 
of the KT line) and continuing to infinitely, positive or 
negative, large values of /14. They map into each other 
with changing the sign of h^. The latter are lines of 
second-order phase transitions with conventional power 
law singularities characterized by rj = 1/4. The exponent 
v diverges at the confluence of these critical lines as 



where y is a parameter that appears in the Hamiltonians 
defined in the range [0, 4/3] for the AT |EJ and [0, 2] for 
the 8V (42). The 0(2) multicritical point identified by 
Kadanoff is the F-model limit of the 8V that is charac- 
terized by y = 3/2, allowed to the 8V but forbidden for 
the AT. 

Another class of models related to the Landau- 
Ginzburg Hamiltonian (II]) is the discrete-cubic N- 
component model |2^{44| It is a short-range inter- 
acted system with A-component spins Si pointing to 
the faces of an A'-dimensional hypercube (face-centered- 
cubic model). The Hamiltonian may be written as: 



n 



^(CSauajSiSj 



(8) 



where is a Potts-like variable that determines which 
component of Si is non-zero and Si an Ising variable that 
determines the sign of that component. For P = the 
Hamiltonian (||) reduces to the 2A-state Potts model, 
for C — to two decoupled A^-state Potts model, and for 
N = 2 to the AT model. The continuous cubic model 
([j]) reduces to the discrete one (|s|) in the limit of strong 
anisotropy (\v\ 3> \u\). In Ref. |46[ it is shown that the 
iteration of RG transformations enforces the continuous 
model, with vq < 0, to have spins pointing only to the 
faces of the hypercube, and so, also for finite anisotropy, 
the two models are equivalent at criticality. The model 
(||) exhibits four competing possible types of critical be- 
havior, related to the Ising model, the A^-and 2A^-state 
Potts models and to a "cubic" fixed point . In Refs. [E_7 48 
it was found that the critical behavior of the discrete 
face-centered cubic model is 0(A^)-like for N < N c = 2, 
AT- like for A" = 2, and characterized by a first-order 
phase transition for N > 2. This result is not surprising 
since the model (||) is related to (|l|) in the region with 
vq < where we expect a first-order phase transition for 
N > N c — 2 and 0(N) behavior for N < 2. 
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In the region with vq > the continuous model (|l|) is 
related to the corner-cubic model in which TV-component 
spins Si point to the corners of an iV-dimensional hyper- 
cube. This model is equivalent to the face-centered one 
for A — 2 (this is the symmetry (||) for the Hamiltonian 
(|l|) and the changing of the sign of hi for (|J)). For other 
values of A the critical behavior of the face-centered- and 
the corner-cubic model is in principle very different. 

This difference is clarified by the study of A-color 
Ashkin- Teller model, first introduced by Grest and 
Widom |4^] , that is equivalent to the corner-cubic model. 
The Hamiltonian of this model reads 
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N 



' N 



(9) 



where sf = ±1 for a — 1, . . . , N are A Ising variables 
and the sum is only over nearest-neighbor (the sign of 
J4 is the opposite of uq in Eq. (Q)). In Ref. it is 
concluded that the order Jj in the RG equations makes 
the decoupled Ising fixed point stable for perturbations 
with J 4 < (u > 0) and unstable for J4 > (u < 
0). In the latter case the system flows away towards 
more negative values of u and finally the transition is first 
order. This has been confirmed jl9| by a Monte Carlo 
simulation for A — 3. This argument was proved exactly 
in the A — > 00 limit [^0|,^l). After, using a mapping of 
the Hamiltonian (O) onto the 0(N) Gross- Neveu model 
p2[ it was argued jp2|-|55|] that the critical behavior of the 
Ising fixed point is affected by logarithmic corrections. 

Finally the three couplings (A, 2) model of Domany 
and Riedel @ (equivalent for N = 3 to the Z{6) 
model pq]) reduces to the A-component discrete cubic 
model, for particular values of the parameters entering 
in the Hamiltonian. The model (3, 2) was studied in 
Refs. [ fl6| , j56| with the Migdal-Kadanoff Renormalization 
Group approach (MKRG). Unfortunately the MKRG is 
expected to give neither precise values of the exponents 
nor the correct nature of the transition, but only a good 
description of the phase diagram. 



C. Order-disorder transition in adsorbed monolayers 

We have already mentioned that the Hamiltonian (|l|) 
describes some order-disorder transitions in adsorbed 
monolayers. The great interest in these models is jus- 
tified since they provide a unique possibility to study ex- 
perimentally a rich variety of two-dimensional systems. 

Several experimental works have confirmed that these 
transitions belong to the universality class of the Ising 
|57-f59f , the three- [|60|]6ll and four-state Potts models 
p2,|6l|, and the XY model with cubic anisotropy [|63|-|66|, 
according to their classification jl|-[|. Monte Carlo sim- 
ulations and other numerical works confirm this scenario 
i7H69|. 



According to the standard classification, the criti- 
cal behavior of the order-disorder transition of diatomic 
molecules on honeycomb lattice in the p(2 x 1) structure 
is described by the A = 3 cubic model dE^l • The ex P er ~ 
imental investigation [[7Q| of Oxygen on Ru(001) shows 
that this transition has critical exponents in agreement 
with the three-state Potts model within about 10%. Sev- 
eral numerical simulations (after the first ones favoring a 
first order phase transition ]7l]j72|| ) confirm this critical 
behavior |73|,[74| . So it may be possible that a fixed point 
of the Heisenberg model with cubic anisotropy has the 
critical behavior of the three-state Potts model. 

We want still to note the possibility for systems having 
exponents close to the four-state Potts model, to belong 
to the universality class of the XY model with cubic 
anisotropy and living near the four-state Potts transition 
characterized by y = 0. 



III. THE FIXED DIMENSION PERTURB ATIVE 
EXPANSION IN TWO DIMENSIONS 

A. Renormalization of the theory 

The fixed-dimension field-theoretical approach fjlf l 
represents an effective procedure in the study of the crit- 
ical properties of systems belonging to the O(N) univer- 
sality class (see, e.g., Ref. |jT5| ). The idea is to extend this 
procedure to models where there are two 4 couplings 
with different symmetry |22] , |l8|| . One performs an expan- 
sion in powers of appropriately defined zero-momentum 
quartic couplings and renormalizes the theory by a set 
of zero-momentum conditions for the (one-particle irre- 
ducible) two-point and four-point correlation functions: 



0(P 4 )] 



(10) 



r ilL(°) = V™ 2 [f + s ^ Sbd + s <*&bc) + 

v 5 ab 5 ac S ad ] . (11) 

They relate the second-moment mass m, and the zero- 
momentum quartic couplings u and v to the correspond- 
ing Hamiltonian parameters r, Uo, and vq: 



u = m uZ u Z, , 



v = m 2 vZ v Z^ 2 . (12) 



In addition, one introduces the function Z t defined by 
the relation 



r^ 2 >(o) = <w?r\ 



(13) 



where r^ 1 ' 2 ) is the (one-particle irreducible) two-point 



function with an insertion of \4> 2 ■ 



From the pertubative expansion of the correlation 
functions r' 2 ', T^ 4 ', and r' 1 ' 2 ) and the above relations, 
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one derives the functions Z<p(u, v), Z u (u, v), Z v (u, v), and 
Z t (u, v) as a double expansion in u and v. 

The hxed points of the theory are given by the common 
zeros of the /3-functions 



P u {u, v) = m 
/3 v (u, v) = m 



du 
dm 

dv 
dm 



u ,v 



(14) 



The stability properties of the fixed points are controlled 
by the eigenvalues u>i of the matrix 



n = 



dj3 u {u,v) d(3 u (u,v) 

du dv 
d/3 v (u,v) d/3 v (u,v) 



(15) 



dv 



dv 



computed at the given fixed point: a fixed point is stable 
if both eigenvalues are positive. The eigenvalues u>i are 
related to the leading scaling corrections, which vanish 
as ~ \t\ Ai where A, = i/u>j. 
One also introduces the functions 



r)<t>(u,v) 
Vt(u,v) 



din Zj, 



9 In to 
dhiZ t 



d In m 



Pu 



= A 



u ,v 



dlnZ^ 

dv 

dlnZ t 
du 



f Pv 

Pv 



dv 
d\nZ t 



du 



Finally, the critical exponents are obtained from 

v = [2-r ]ck (u*,v*) + r) t (u*,v*)]- 1 , 
7 = v{2 - rj), 



(16) 
(17) 



(18) 
(19) 
(20) 



where (u* , v*) is the position of the stable fixed point. 



B. The four loop series 

In this section we present the perturbative expansion 
of the RG functions ([l4|), (|l6|), and ( |l7| ) up to four loops. 
The diagrams contributing to the two-point and four- 
point functions are reported in Ref. J7(|. We do not cal- 
culate the integrals associated to each diagram, but we 
use the numerical results compiled in Ref. |76| . Summing 
all contributions with the right symmetry and group fac- 
tors (see Ref. [§) we obtain all RG functions. 

The results are written in terms of the rescaled cou- 
plings 



8tt 



Rn u, 



8tt 



v, 



(21) 



where Rn = 9/(8 + TV). We adopt this rescaling to have 
finite fixed point values in the limit N — > oo. 
The resulting series are 



/3 fi = -u + u 2 + -uv + u E &y ( 22 ) 

i+j>2 



ft, = -v + v 2 + ^7^™ + v E (23) 



V<t> 



E 

i+j>2 



111 



2{2 + N) r 
' (8 + N) '' 



E 

i+j>2 



where 



3 

loTT 



^2nP u ^ 



Pv 



3 
16^ 



(24) 



(25) 



(26) 



The coefficients b. 



(«) u(v) W 



b)V , e)V , and ef^ are reported in 



the Tables Jj, y, Q and |v|. Note that due to the 
rescaling (|2f), the matrix element of Q are two times 
the derivative of (3 with respect to u and v. 

We have verified the exactness of our series by the fol- 
lowing relations: 

(i) f3u(u,0), 770(72,0) and 77* (zZ, 0) reproduce the corre- 
sponding functions of the 0(iV)-symmetric model {77].[78| : 

(ii) j3v(0,v), ?70(O, v) and r]t(0,v) reproduce the corre- 
sponding functions of the Ising-like (N — 1) </> 4 theory; 

(iii) The following relations hold for N = 1: 



j3 a (u, x — u) + j3y(u, x — u) = /3 B (0, x), 
ri<p(u,x- u) = 770(0, x), 
rj t (u,x- u) = ri t (0,x). 



(27) 



(iv) For N = 2, using the symmetry (^|) and ([|), and 
taking into account the rescalings (|2l]), one can easily 
obtain the identities 

5 - 

Pu(u+ |tJ,-7j) + -Pv[u+jV,-v) = Pu(u,v), (28) 
(3y(u + |tj, -tj) = -j3v(u,v), 

?70(M-|- §7J,-7j) = 7/0(71,7;), 



\v, -v) = rj t (u,v). 



These relations are exactly satisfied by our four-loop se- 
ries. Note that, since the Ising fixed point is (0,g£)> an d 
g*j is known with very high precision from Ref. |79j 



g*j = 1.7543637(25), 



(29) 



the above symmetry gives us the location of the cubic 
fixed point: (| 5 J,- 5 |). 

(v) In the large-TV limit the critical exponents of the 
cubic fixed point are related to those of the Ising model: 
77 = rji and v = fj . One can easily see that, for N — > 00, 
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2fl -(10.3350 + 47. 675TV) 

1,1 -8.39029 

0,2 -0.21608 

37) 524.377 + 149.152TV + 5.00028TV 2 

2.1 144.813 + 7.277557V 

1.2 10.0109 + 0.05832787V 

0,3 0.231566 

4~0 -(7591.108 + 2611.157V + 179.697TV 2 + 0.0888427TV J ) 

3.1 -(2872.09 + 291.2557V - 0.1268147V 2 ) 

2.2 -(330.599 + 5.970867V) 

1.3 -(16.0559 + 0.05789557V) 
4,0 -0.311695 



2,0 


6.751258 (TV + 2) 


1,1 


4.500839 


0,2 


0.250047 


3,0 


-(96.7105 + 65.16867V + 8.406687V 2 ) 


2,1 


-(48.3553 + 8.406687V) 


1,2 


-(6.19023 + 0.1166567V) 


0,3 


-0.233588 


4,0 


1135.05 + 844.57V + 139.6567V 2 + 0.5833777V J 


3,1 


756.697 + 184.6527V + 0.7778367V 2 


2.2 


149.467584 + 7.5534597V 


1,3 


11.515397 + 0.1157917V 


0,4 


0.323089 



TABLE I. The coefficients cf. Eq. (|22|) 



TABLE IV. The coefficients ef}, cf. Eq 



(TV + 8) i by! 



TABLE II. The coefficients cf. Eq 



1 ( 1 
0,2 

3,0 
2,1 
1,2 
0,3 
4,0 
3,1 
2,2 
1,3 
0,4 



(TV + 8)'fr<; 



-(92.6834 + 5.83417TV) 

-17.392 

-0.716174 



1228.63 + 118.504TV - 1.83156TV 2 

358.882 + 2.84758TV 

31.4235 

0.930766 

-(20723.1 + 2692. OTV + 25.48547V 2 + 0.8246557V J ) 

-(8273.28 + 233.787V - 0.5747577V 2 ) 

-(1134.8 + 1.914027V) 

-68.4022 

-1.58239 



TABLE III. The coefficients e\j', cf. Eq.i 



h J 

1,1 
0,2 
3,0 
2,1 
1,2 
0,3 
4,0 
3,1 
2,2 
1,3 
0,4 



(TV + 8)'e 



0.917086 (TV + 2) 

0.611391 

0.033966 



■i 



-0.054609 (TV - 
-0.054609 (TV - 
-0.054609 
-0.002023 



2)(TV + i 



(TV + 2)(26.7676 + 4.24178TV - 0.092684TV 2 

27.6901 + 5.6557TV - 0.123579TV 2 

5.404243 + 0.132800TV 

0.410151 

0.011393 



7?0 {u,v) — f]j(v), where rjj(v) is the perturbative se- 
ries that determines the exponent r\ of the Ising model. 
Therefore, the first relation is trivially true. On the other 
hand, the second relation v = vi is not identically satis- 
fied by the series, and is verified only at the critical point 
& 

(vi) For N = the series reproduce the results of Rcf. 



C. Resummations of the series 

The field theoretic perturbative expansion generates 
asymptotic series that must be resummed to extract the 
physical information about the critical behavior of the 
real systems. 

Exploiting the property that these series are Borel 



summable for <ji 4 theories in two and three dimensions 
pfj| , one can resum these perturbative expressions con- 
sidering the Borel transform combined with a method for 
its analytical extension. In the case of the O(N) symmet- 
ric model with only a coupling g, all perturbative series 
are of the form F{g) ~ fkg k ■ Exploiting the knowl- 
edge of the large order behavior of the coefficients fk (Ref. 



fk ~ k\ (~a) k k b [l + Oik- 1 )] 



with 



(30) 



(a large order behavior related to the singularity of the 
Borel transform closest to the origin) one can perform the 
following mapping pV( 



y{g) 



V 1 -ff/fffc- 1 
y/1 - g/g b + 1 



(31) 



to extend the Borel transform of F(g) to all positive val- 
ues of g . The singularity g\, depends only on the consid- 
ered model and can be obtained from a steepest-descent 
calculation in which the relevant saddle point is a finite- 
energy solution (instanton) of the classical field equations 
with negative coupling J8^,Q. Instead the coefficient b 
depends on which Green's function is considered. 

Note that the function F(g) can be Borel summable 
only if there are no singularities of the Borel transform 
on the positive real axis. 

This resummation procedure has worked successfully 
for the O(N) symmetric theory, for which accurate esti- 
mates for the critical exponents and other physical quan- 
tities have been obtained ]8l]j8^,^5| . 

For this reason we want to extend the resummation 
procedure cited above to multicoupling models, as it has 
been done for the three-dimensional cubic model |L8|j78| 
and for the frustrated system with non-collinear order in 
two dimensions Considering a double expansion in 
u and v at fixed z = v/u, and studying the large or- 
der behavior (following the same procedure used in Refs. 
|l8|,^6|) of the new expansion in powers of u to calculate 
the singularity of the Borel transform closest to the origin 
Ub we have 



G 



— =-a{R N + z) for < z, (32) 

Ub 



1 

Ub 



— = — a (Rn + jjZ) for > z > — 



2A^i?A, 
7V + 1 ' 



where a = 0.238659217 .... 

Note that the series in powers of u keeping z fixed is not 
Borel summable for u > and z < —Rn- This fact will 
not be a real limitation for us, since we will only consider 
values of z such that Ub < 0. It should be noted that 
these results do not apply to the case N = 0. Indeed, in 
this case, additional singularities in the Borel transform 
are expected J87j . 

The exponent 6 in Eq. (^5|) is related to the number 
of symmetries broken by the classical solution |Q. It 
depends on the quantity considered. In the cubic model, 
for v 7^ 0, we have 6=2 for the function 77^, and 6 = 3 for 
the /3-functions and ?y t . For v = 0, we recover the results 
of the 0(A r )-symmetric model, that is 6 = (3 + N)/2 for 
77^, and 6 = (5 + N)/2 for the /^-function and r) t Q. 

Good estimates of the critical exponent could be ob- 
tained also using a Pade Borel analysis of the series, as 
shown in the case of the two-dimensional O(N) models 
and the random Ising [jL5|. 

An important issue in the fixed dimension approach to 
critical phenomena (and in general of all the field theoret- 
ical methods) concerns the analytic properties of the (3- 
functions. As shown in Ref. [^9j for the 0(N) model, the 
presence of confluent singularities in the zero of the per- 
turbative (3 function causes a slow convergence of the re- 
summation of the perturbative series to the correct fixed 
point value. The O(N) two-dimensional field-theory esti- 
mates of physical quantities J8|,|7^] are less accurate than 
the three-dimensional ones, due to the stronger nonana- 
lyticities at the fixed point |^]-|92| . In Ref. Q it is shown 
that the nonanalytic terms may cause large imprecisions 
in the estimate of the exponent related to the leading cor- 
rection to the scaling u>; instead the result for the fixed 
point value is a rather good approximation of the correct 
one (if one compare the field theoretical results for the 
four point renormalized coupling in the g-expansion |7§| ] 
and in the e-expansion p2] , p3[ w ith the estimates of other 
non-perturbative methods ]79| , p4| , p5| , one finds that the 
systematic error is always less then 10%). We think that 
this scenario holds also for the cubic models. 



IV. FOUR-LOOP EXPANSION ANALYSIS 

A. The analysis method 

In order to study the critical properties of the continu- 
ous anisotropic cubic model (Q), we use two different re- 
summation procedures: the Pade-Borel method and the 
conformal mapping of the Borel-transformed series. 



Explicitly, let us consider an Hoop series in u and v of 
the form 



I l-k 



[z u 



(33) 



k=0 h=0 



k=0 



where R is one of the RG functions and z is the ratio v/u 
that we will consider always fixed to same value. In this 
manner we have an asymptotic series of only one variable 
u, depending on the additional parameter z that have to 
run from to 00 in order to reproduce all the quadrant 
with u, v > 0. In order to use a variable defined in a 
finite range we set 



7T 

tan —x, 



(34) 



with x € [0, 1]. In the follow we fix x — k/N max , where 
k is an integer number running from to N max (we will 

fix iV moa! = 25). 

In the Padc-Borcl method we consider the Pade 
approximants [L/M] of the Borel-Leroy transform of 
R(u, z). Explicitly the approximants of R arc 



E(R)(L,M,b;u, z) 



dtt b e~ 



■ N L (ut,z) 
D M (ut, z) 



(35) 



where Nl and Dm are two polynomials of degree L and 
M respectively, with L + M < /, that are determined by 
the condition that the expansion of E(R) (L, M, 6; u, z) in 
powers of u gives R(u, z) to order L + M . 

In this manner we have several approximants of the 
function R(u, v) with varying the three parameters 6, L, 
and M. As usual in this case |Q the best estimates 
of the resummed function are given by the diagonal and 
near-diagonal approximants (i.e. for our four-loop series 
we have three reasonable choices [2,2], [3, 1], and [2, 1]). 
Then we search the value of 6 (called b opt ) integer or 
half-integer minimizing the differences between the con- 
sidered approximants; a reasonable estimate will be the 
mean value of E(b opt , L, M) on all the values of L and M 
considered. Then we could take as error bar the devia- 
tions from the mean value of all the approximants with 
b pt — 1 < b < b opt + 1 (always considering integer and 
half-integer values of 6). 

In the analysis using the conformal-mapping method 
we essentially follow the procedure used in Ref. . We 
exploit the knowledge of the value of the singularity of the 
Borel transform closest to the origin (a value given in the 
previous section), and we generate a set of approximants 
to our asymptotic series, varying the two parameters a 
and 6 appearing in 



E(R) p (a, 6; u, z) = } B k (a 1 6; z) 



k=0 



b y(ut;z) k 
1 [1 -!,(«<;*)]«' W 
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TABLE V. Half of the exponent lo 2 at the O(N) fixed 
point. CM is the value obtained using conformal mapping 
technique and PB the one using Pade-Borel. 



N 


u O(N) 


CM 


PB 


2 


1.80(3) 


0.03(3) 


0.06(4) 


3 


1.75(2) 


-0.08(3) 


-0.07(3) 


4 


1.70(2) 


-0.18(4) 


-0.17(5) 


8 


1.52(1) 


-0.45(5) 


-0.44(6) 



where 



y{x;z) = 



x/ub{z) — 1 
V'l - x/u b (z) + 1 



(37) 



The coefficients Bk are determined by the condition that 
the expansion of E(R) p (a,b;u, z) in powers of u and v 
gives R[u, v) to order I. 

The procedure to choose the range of the parameters 
a and b used to find good estimates and reasonable error 



bars of the various quantities is the same of Rcf. 18 



B. Stability properties of the O(N) and the Ising 
fixed points 

First of all, we analyze the stability properties of the 
0(-/V)-symmetric fixed point. Since 



¥ (M) = ' 



(38) 



the eigenvalues are simply 

^1 = 2— (u ,0), 
ou 



lo 2 = 2^(u*,0), (39) 
ov 



where u* is the fixed-point value of the O(N) vector 
model of which accurate estimates are available from the 
five-loop analysis of Ref. |7q| . The exponent u>i is the 
usual exponent of the 0(A r )-symmetric theory |7a] that 
gives the first correction to the scaling. For all TV > 2 it 
is known that uj\ assumes the constant value 2 (see Ref. 
|39"f and references therein) . Instead u>2 is the eigenvalue 
that determines the stability of the fixed point respect to 
an anisotropic cubic perturbation. 

In Table [v| we report the results for ui2 for several val- 
ues of N. It is quite evident that the O(N) fixed point 
is unstable for N > 3 (we note that for small N these 
values are very close to zero). In the limit N — > oo it 
holds u>2 — —2. For N = 2 our result is compatible with 
the presence of a marginal operator, i.e., lu2 = 0. 

Then we focus our attention on the stability properties 
of the Ising fixed point. Also in this case the eigenvalues 
are simply: 



2-^(0, « ), 



LU 2 



2-—(Q,v ), 
ou 



(40) 




0.5 1 1.5 2 

FIG. 2. Zeros of the /3-functions for iV = 2 in the (u, v) 
plane using the conformal mapping method. The straight 
line correspond to zeros of /3 a (u, v) and the dashed one to the 
zeros of (3v(u, v). 



where v* is the fixed-point value of the Ising model (|2 
The exponent cJi gives the correction to the scaling of the 
Ising model and it should be equal to 7/4 Q, while u>2 
is the eigenvalue determining the stability of this fixed 
point. We find that the series u>2(v) is independent from 
N (we already know that its fixed point value must be 
equal to ai for all N): 



uj{{v) 



-1 + -V- 0.21617J 2 + 0.231577J 3 



0.311697J 4 . 
(41) 



The fixed point value of this exponent is w|/2 = 
—0.10(5), using the conformal mapping method, and 
—0.08(5), using the Pade-Borel analysis. These values 
are quite close to —aj — 0; we attribute this small dis- 
crepancy to the presence of nonanalyticities at the zero 
of the (3y. To support this thesis we note that the value 
of W2/2 is approximately the same of the resummation of 
the exponent —aj = —0.081 found in Ref. |7q|. 



C. Evaluation of non-trivial fixed points 

Since the Ising fixed point has always a marginal op- 
erator and the O(N) is unstable for N > 3, we search 
for the presence of non-trivial fixed points in the region 
U,V > 0. 

In Figs. H and we show the zeros of the /? functions 
for N=2 using the conformal mapping technique and the 
Pade-Borel analysis respectively. The curves of zeros of 
/3u and /3& seem parallel and do not have any intersec- 
tion (in these figures, and in all the others, the error bar 
are not shown to make them more readable). The dis- 
tance between the two lines is always less than the error 
bar (that is of the order of 0.1 for each line). So our re- 
sult is compatible with the presence of the line of fixed 



8 



points conjectured by Kadanoff joining the decou- 
pled Ising and the 0{2) fixed points. In the Kadanoff 
identification this line of fixed points is the dual line of 
the 8V model in the region that goes from the F-model 
limit (0(2) multicritical point with y — 3/2) to the de- 
coupled Ising point (y = 1). This line of fixed points is 
identifiable with the one of Ref. (3{| with hi < 0. 

For the symmetry (||) the fixed-point line continues 
for negative v and reaches the "cubic" fixed point (that 
in this case is a standard Ising) located at (5/3 g'j, —g})- 
This second line is the one of Ref. Q with > that is 
isomorphic to the other one with /14 < 0. We do not know 
if the line continues in the region with negative u, that 
may be a region of second order phase transition, having 
the critical exponent of the 8V model with y < 1 (only 
for v > —u, from the stability condition). We tried to 
perform the resummation also for these value of u, v but 
we obtained strongly oscillating results. 

For the other values of N we find one fixed point at the 
intersection of the two curves of zeros of the /3-functions, 
that is the usual cubic fixed point, analytic continuation 
of the one found in the e-expansion. For example the 
results of the conformal mapping analysis for N = 3 and 
of the Pade-Borel for N — 4 are shown in Figs. || and || 
respectively. All our results about the cubic fixed points 
are summarize in Table VI. We note that, within the er- 



ror bar, there is full agreement between the two methods 
& 

In the limit N — ► 00 the series (3u(u, v) and (3^(u, v) 
simplify to 



Aj(u,u)^ /% s ing( u )' 
Ar(u, u) -> u[Pi(v) - uP 2 {v)] 7 



(42) 
(43) 



with 




0.5 1 1.5 2 

FIG. 4. Zeros of the /3-functions for N = 3 in the (u, v) 
plane using the conformal mapping method. The straight 
line correspond to zeros of /3 a (u, v) and the dashed one to the 
zeros of (3v(u, v). 



2 

■l + -v 



0.21608?T + 0.23157F 



—3 



0.311681T 4 



and 



P 2 {v) = l + 0.058328w 2 



0.0578961J 3 . 



(44) 



(45) 



The zero of fiv(u, v) is constant with varying u at the 
value v* = Qj (|2^) , and so the u coordinate of the cubic 
fixed point is simply: 



u JV-*c 



Pi(g*i) 



0.09(4). 



(46) 



We note that the small value of u* for TV — > 00 does not 
exclude that the Ising and the cubic fixed points coalesce 
in this limit. 




0.5 1 1.5 2 

FIG. 3. Zeros of the /3-functions for iV = 2 in the (w, v) 
plane using the Pade-Borel analysis. The straight line corre- 
spond to zeros of /3u(u, v) and the dashed one to the zeros of 

f3v{u,v). 




0.5 1 1.5 2 

FIG. 5. Zeros of the /3-functions for iV = 4 in the (u, v) 
plane using the Pade-Borel analysis. The straight line corre- 
spond to zeros of 0^(u, v) and the dashed one to the zeros of 
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TABLE VI. Cubic fixed points for N > 3. 



N 


CM 


PB 




(«*,«*) 


(«*,«*) 


8 


(0.24(8), 1.72(10)) 


(0.15(9), 1.72(12)) 


4 


(0.54(10), 1.43(8)) 


(0.38(11), 1.51(9)) 


3 


(0.83(12), 1.12(9)) 


(0.71(12), 1.17(9)) 



TABLE VII. The half of eigenvalues of the Q matrix for 
TV > 3 at the cubic fixed point. 



We want to add that in Ref. Q, using an argument 
based on conformal field theory, it is conjectured that the 
cubic fixed point and the Ising one merge for arbitrary 
N. Our analysis seems to contradict this statement. 



D. Stability properties of the cubic fixed point 

For N = 2 we check that one eigenvalue of the stabil- 
ity matrix Q vanish on all the line of fixed points. From 
Figs. H it is evident that this line is almost straight, so 
for our calculation we evaluate the f2 matrix on the line 
joining the Ising fixed point (at (0, <?J)) and the 0(2) (at 
(1.80(3), 0) from Ref. [78)). Within the precision of our 
calculation we confirm that on this line one marginal op- 
erator exists (see Fig. ||, the uncertainty of w 2 /2 is con- 
stantly of the order of 0.05). 

The results for the eigenvalues of the f2 matrix at the 



cubic fixed point for N > 3 are summarize in Table VII. 
The eigenvalue L02 is positive for all considered values of 
N, but it is very close to zero. Being the Ising fixed 
point stable against perturbations with positive u, it is 
not possible that the cubic fixed point is stable for "left" 
perturbations. We believe that our results are compatible 
with the presence of a marginal operator also at the cubic 
fixed point and that the very small discrepancy of 0J2 with 
respect to zero is due to nonanalytic terms. 

We conjecture that the RG flow in the coupling plane 
U, v for N > 3 is the one sketched in Fig. [?]. For 



w/2 
. 8 r 



0.6 



0.4 



0.2 




N 


CM 


PB 


pseudo-e 




(wi/2,w 2 /2) 


(wi/2,w 2 /2) 


(wi/2,w 2 /2) 


00 




(0.64(10), 0.11(8)) 


(0.65(3), 0.07(4)) 


8 




(0.70(10), 0.07(6)) 


(0.72(3), 0.08(4)) 


4 


(0.65(15), 0.11(6)) 


(0.70(15), 0.06(6)) 


(0.79(4), 0.06(4)) 


3 


(0.70(7), 0.09(6)) 


(0.70(10), 0.04(5)) 


(0.75(6), 0.05(5)) 



U, v < we expect a first order phase transition. For 
< v/u < v* /u* the transition is second order and its 
critical behavior is characterized by the cubic fixed point. 
For v/u > v* ju* the transition is governed by the Ising 
fixed point, i.e., with the exponents of the Ising model 
more logarithmic corrections. 

To support this scenario we also consider the pseudo-e 
expansion using substantially the same procedure of Ref. 
fl27| . By using this trick we have a smallest error bar since 
we avoid the uncertainty of the fixed point (for a detailed 
discussion see Ref. ) . The res ults are equivalent to the 
ones presented above (see Table VII ). 

The presence of a marginal operator both at Ising and 
at cubic fixed points means that we should expect slowly 
decaying crossover effects that could produce systematic 
errors in Monte Carlo simulations and in the experimen- 
tal investigations. 



E. Critical exponents 

The direct evaluation of the critical exponents for the 
two-dimensional O(N) model leads to erroneous values 
because of the strong effect of nonanalytic terms. In fact 



FIG. 6. Values of ui\/2 (dashed line) and 0*2/2 (straight 
line) with varying the parameter x (B4J). 




FIG. 7. Conjectured renormalization-group flow in the 
coupling plane (u, v) for N > 2. 
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from the analysis of the four- and five-loop series of the 
N = 1 model is found r) ~ 0.131 @ and r) = 0.13(7) 
plf instead of the exactly known r\ = 1/4. The value of 
v for the Ising universality class is quite good (y ~ 0.96 
|pl|,[r8[ ). An our unpublished analysis of the O(N) se- 
ries (already cited in Ref. |8(|) shows that for higher val- 
ues of N the effect of nonanalyticities is very dangerous. 
In fact, studying four-loop series, we find r\ = 0.11(6) 
both for the XY and Heisenberg model instead of 1/4 
and 0. In the same way for rj — r\ t we find 1.18(5) and 
1.36(4), for N = 2, 3, instead of 2. This latter systematic 
error brings to a finite value of v. The above standard 
analysis applied to the XY model with cubic anisotropy 
gives 77 ~ 0.11 along the fixed-point line, and values rang- 
ing from 1.18 to 0.97 for 77 — r/t. 

In order to reduce the effect of nonanalytic terms in the 
estimates of the critical exponents of the cubic model we 
adopt a new strategy. We use a constrained analysis on 
the two variables series 77 and 77 — r\ t fixing the values 
assumed at the fixed point on the axes u = and v = 0. 
Explicitly, let us consider a generic /-loop series in u and 
v without constant term 



( l-k 



R(u,v) = J2J2 Rhkllh v k = J2 Rk & 



(47) 



fc=0 /i=0 



k=l 



of which we know the values R(u* , 0) = a and i?(0, v*) — 
b. We can rewrite the previous function in the form 



R(u,v) = R(u,0) + R(0,v) + uvA R (u,v) 



(48) 



where uvAr(u, v) is the difference between the orig- 
inal function R(u,v) and the value that it as- 
sumes on the axes. We resum the three functions 
R(u,0) , i?(0,?7), and Ar(u,v), appearing in (fl8|), in an 
independent way. In the case of R(u, 0) and R(0,v) we 
adopt the standard method to resum one variable func- 
tions with some constraints, as done for the e-expansion 
in Refs. A full description of this method may 

be found in Ref. @. 

First of all we consider the peculiar case N = 2. We 
evaluate the critical exponents on the straight line join- 
ing the Ising and the 0(2) fixed points (as in the case of 
the f2 matrix). The results for the exponents 77 and f]—r\t 
are presented in the Figs. || and ^ respectively. Within 
the precision of our calculation we find an exponent 77 as- 
suming the constant value 1/4 along the line as predicted 
in Ref. ]35|] (the uncertainty is about 0.05 for x ~ 0.5 and 
decreases near the borders) . The exponent 7/ — 77 f inter- 
polates in a monotonic and continuous way from 1 at the 
Ising fixed point to 2 at the 0(2) fixed point and so it 
is compatible with the Kadanoff's conjecture identifying 
this line with the dual line of the 8V model with the 
continuous varying exponent given by (Q). 

The parameter y in equation (]?]) belongs to the range 
[1,3/2] respectively from Ising to 0(2), instead x ranges 
between and 1 from 0(2) to Ising. So we try to identify 



0.4 

n 

0.35 
0.3 

0.25 
0.2 

0.15 
0.1 

0.05 



0.2 0.4 0.6 0.8 1 

FIG. 8. Values of 77 with varying the parameter x ([S]). 
The straight line represents CM results and the dashed the 
PB ones. 



y = 



(49) 



Using this assumption the exponent 77 — 77* could be writ- 
ten in terms of x as 



V~Vt 



l + x 



(50) 



This curve is shown in Fig. ^| together with the numeri- 
cal values obtained from the resummation of the pertur- 
bative series: all the curves are very close. This result 
strongly supports both Kadanoff's conjecture about the 
nature of the line of fixed points and the relation ( |49| ) 
between the parameter y of the 8V-model and x of the 
XY model with cubic anisotropy. 

The exponents at the cubic fixed point obtained using 
this constrained analysis for some N > 3 are reported in 
Table VIII. We find an 77 exponent varying very slowly 



with TV. Our estimates are also compatible with an 77 
constant and equal to 1/4 for all N. The estimates of 




0.2 0.4 0.6 0.8 X 1 

FIG. 9. Values of 77 — r\ t (straight line CM and pointed line 
PB) with varying the parameter x (^). The dashed line is 
our conjecture (ptl). 
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rj — rjt are instead affected by a big error bar, mainly due 
to the uncertainty of the cubic fixed points. We could 
reduce the last error using the pseudo-e expansion, but 
in this way we will find again the systematic error of 
nonanalyticiti es tha t we cannot control. 

From Table VIII it is clear that when N increases the 



exponents get closer to the ones of the Ising model, that 
is the limit for N — > oo, cf. (|J). In this limit the series 
reproduces order by order the one of the Ising model. 
For the difference rj — rj t the equality of the exponent is 
expected only at the fixed point. In fact it holds 



N 



lim rjtiu, v) = —2u + rj((v) — uv 2 A(v) (51) 



where r\\ (v) is the series of the Ising model, and A(v) a 
function that at four-loop reads A{v) — 0.1166 — 0.1157?;. 
At the cubic fixed point we obtain 



nt - vi = o.i3(io) 



(52) 



that is compatible with zero. 

For N = 3 the result for r\ — r\ t via scaling law leads 
to v = 1.3(3), so there are no fixed points (stable or 
unstable) with the three-state Potts exponent v = 5/6. 
In Sec. [I C we have stressed that experimental investi- 
and Monte Carlo simulations [|73|J74| in the 



gations 

adsorption of diatomic molecules on honeycomb lattice 
in the p(2 x 1) structure display the critical behavior of 
three-state Potts, although, according to the standard 
classification § , this system should be described by the 
TV = 3 cubic model. It may be possible that the system 
undergoes a weak first-order phase transition in the face- 
cubic region with v < 0, with effective exponents close to 
the ones of the three-state Potts model which appears in 
the discrete face-centered model. This point needs fur- 
ther studies for a full understanding of the phenomenon. 



V. CONCLUSIONS 

In the present paper we have studied the critical behav- 
ior of TV-component spin models with cubic anisotropy by 
applying the field-theoretic renormalization-group tech- 
nique directly in two-dimensions. 

We have firstly focused our attention to the XY 
model (N = 2) with cubic anisotropy. We found that 
this model has a line of fixed points joining the decoupled 
Ising and the 0(2). Along this line the critical exponent 

TABLE VIII. Critical exponents for N > 3 at the cubic 
fixed point. 



N 


rj [CM.PB] 


V - rjt [CM,PB] 


8 


[0.26(2),0.25(2)] 


[1.20(15), 1.15(15)] 


4 


[0.23(4),0,20(6)j 


[1.22(15), 1.18(20)] 


3 


[0.17(8),0.15(9)j 


[1.39(20), 1.23(20)] 



■q assumes the constant value 1/4. The exponent v runs 
from 1, at Ising, to oo, at 0(2\ according to the standard 
expression of the 8V model (0). This argument is based 
on a mapping of the y parameter of the 8V model onto 
the x parameter of the Hamiltonian ([!]): 



2 v* 
x = 3 — 2y = — arctan — . 

7T u 



(53) 



The predicted expression for v (Eq. (pCf)) is in very 
good agreement with the resummation of perturbative 
series (see Fig. [)]). In this manner we relate the mea- 
sured exponents to the strength of the anisotropy (the 
parameter x) |)7| . The XY model with cubic anisotropy 
for u < could display the critical behavior of the Askin- 
Teller model with y < 1, i.e., the region on the dual line 
from the decoupled Ising fixed point (y — 1) to the four- 
state Potts (y — 0). We tried to check if the line of fixed 
points continues for negative value of m, but the results 
of the resummation for these values of the renormalized 
couplings are strongly oscillating. 

For all other values of N we found one fixed point in the 
region with u,v > 0, that is the usual cubic fixed point, 
analytic continuation of one found in the e-expansion. 
The fixed point is marginally stable and unstable for 
perturbations with v/u < v* /u* and v/u > v* /u* re- 
spectively. The conjectured RG flow diagram is sketched 
in Figure 0. The estimates of this fixed point for sev- 
eral values of N are reported in Table VI and the critical 



behavior is characterized by the exponents of Table VIII 

In particular for the physical relevant case of N 
the cubic fixed point we found 



3 at 



77 = 0.17(8) , 



v = 1.3(3). 



(54) 



The value found for v is different from the one of the 
three-state Potts model. This fact deserves further stud- 
ies about the order-disorder transitions of some adsorbed 
monolayers that should belong to the universality class of 
the Heisenberg cubic model, but that experimentally dis- 
play the exponents of the three-state Potts model. One 
possible scenario is that the system undergoes a very 
weak first-order phase transition with effective exponents 
close to the ones of three-state Potts model. 
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